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Abstract 

in ; 

q \ The dip is a feature in the diffuse spectrum of UHE protons in energy range 1 x 10 18 — 4 x 

\ 10 19 eV, which is caused by electron-positron pair production on CMB photons. Calculated for 

^| power-law generation spectrum with index 7 9 = 2.7, the shape of the dip is confirmed with high 

C/} \ accuracy by data of Akeno - AGASA, HiRes, Yakutsk and Fly's Eye detectors. The predicted shape 

OO ; of the dip is robust: it is valid for the rectilinear and diffusive propagation, for different discretenesses 
in the source distribution, for local source overdensity and deficit, for source inhomogeneities on scale 

| I < 100 Mpc etc. Below the characteristic energy E c « 1 X 10 18 eV the spectrum of the dip flattens 

Q\ ' for both diffusive and rectilinear propagation, and more steep galactic spectrum becomes dominant 

q \ at E < E c . The energy of transition E tr < E c approximately coincides with the position of the second 

Q\ ; knee E 2 k n observed in the cosmic ray spectrum. The critical energy E c is determined by the energy 

^ ; E eq = 2.3 x 10 18 eV, where adiabatic and pair-production energy losses are equal. Thus, position of 

O ■ the second knee is explained by proton energy losses on CMB photons. 
43 ' 

^p" 1 , 1. Introduction 

o ■ 

^ ' The nature of signal carriers of UHECR is not yet established. The most natural primary 

f3 particles are extragalactic protons. Due to interaction with the CMB radiation the UHE protons 

.^J ' from extragalactic sources are predicted to have a sharp steepening of energy spectrum, so called 

X GZK cutoff p. 

There are two other signatures of extragalactic protons in the spectrum: dip and bump [2] - 
[3] . The dip is produced due to p + 7cmb — > p + e + + e~ interaction. The bump is produced by pile- up 
protons which loose energy in the GZK cutoff. As was demonstrated in jH], the bump is clearly seen 
in the calculations for a single source at large redshift z, but it practically disappears in the diffuse 
spectrum, because individual peaks are located at different energies. 

As it will be discussed in this paper, the dip is a reliable feature in the UHE proton spectrum 
(see also [Hj - (Oj)- Being relatively faint feature, it is however clearly seen in the spectra observed by 
AGASA, Fly's Eye, HiRes and Yakutsk arrays (see jTHj - [15! f° r ^ ne data)- We argue here that it 
can be considered as the confirmed signature of interaction of extragalactic UHE protons with CMB. 

The dip has low-energy and high-energy flattenings (see FigEJ). The high-energy flattening 
at E a ~ 1 x 10 19 eV reproduces the ankle, the feature well seen in the observational data. The 
low-energy flattening at E c ~ 1 x 10 18 eV, seen for both rectilinear [Oj and diffusive [TT)| IT7| ITH] 
propagation, is a natural place for transition from extragalactic to galactic component of cosmic 
rays. Somewhere below E c the more steep galactic component must appear. Thus E c describes the 
beginning of transition (if one follows it from high energy side) or the energy where this transition 
completes (if one moves from low to high energy). We shall demonstrate here (see also [T^j) that E c 
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is connected with E C(l = 2.3 x 10 18 eV, where pair-production and adiabatic energy losses are equal. 
The visible transition from galactic to extragalactic cosmic rays occurs at E tr < E c , and this energy 
coincides well with position of the second knee (Akeno - 6 x 10 17 eV, Fly's Eye - 4 x 10 17 eV, HiRes - 
7 x 10 17 eV and Yakutsk - 8 x 10 17 eV). 

The transition from galactic to extragalactic cosmic rays at E c ~ 1 x 10 18 implies the dominance 
of the proton composition of the observed cosmic rays at E > 1 x 10 18 eV. While HiRes [01] HiRes- 
MIA j2D] and Yakutsk j2Ij data support this prediction and Haverah-Park ^3] data do not contradict 
it at E > (1 — 2) x 10 18 eV, the Akeno [22] and Fly's Eye ^2] data favour the mixed composition, 
dominated by heavy nuclei (for a review see [THj and |23j). 

Below we shall analyze the features in UHE proton spectrum using basically two assumptions: 
the uniform distribution of the sources in the universe and the power-law generation spectrum. We 
do not consider the possible speculations, such as cosmological evolution of sources. Only in Section 
^]we shall turn to model-dependent consideration. 



2. The dip 

The analysis of the dip is convenient to perform in terms of the modification factor. 
Modification factor is defined as a ratio of the spectrum J P {E) with all energy losses taken into 
account, and unmodified spectrum Jp nm , where only adiabatic energy losses (red shift) are included. 



T)(E) 



J P (E) 



Tunm 



Modification factor is less model-dependent quantity than the spectrum. In particular, it depends 




Fig. 1. Modification factor for the power-law generation spectra with 7 9 in the range 2.0 - 2.7. Curve 
r] = 1 corresponds to adiabatic energy losses, curves r] ee - to adiabatic and pair production energy losses 
and curves r\tot - to total energy losses. 



weakly on j g , because both numerator and denominator in Eq. (Q) include E~ la . In this paper we 
shall consider the non-evolutionary case m = 0. In Fig El the modification factors are shown as a 
function of energy for two spectrum indices 7 9 = 2.0 and 7 9 = 2.7. They do not differ much. 
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Comparison of the predicted dip {j] ee curve) with observational data are shown in Fig ED for 
7 9 = 2.7. One can observe the excellent agreement with AGASA, HiRes and Yakutsk data (and also 
with Fly's Eye data not presented here), while the Auger spectrum, at this preliminary stage, does 
not contradict the dip. 

The systematic errors in energy determination of existing detectors exceed 20%. The dip can 
be used for energy calibration of detectors. Assuming the energy-independent systematic errors we 
shifted the energies of AGASA and HiRes by factor AA g = 0.9 and factor Am = 1.2, respectively, to 
reach the best fit of the dip. In Fig ED the spectra of Akeno- AGASA and HiRes are shown before and 
after this energy calibration. One can see the good agreement of the data at E < 1 X 10 20 eV and 
their consistency at E > 1 x 10 20 eV. This result should be considered together with calculations |2*H] . 
where it was demonstrated that 11 superGZK AGASA events can be simulated by the spectrum with 
GZK cutoff in case of 30% error in energy determination. We may tentatively conclude that existing 
discrepancy between AGASA, HiRes and Auger spectra at all energies may have statistical origin. 
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Fig. 2. Predicted dip in comparison with AGASA, HiRes, Yakutsk and Auger data. 



2.1 Robustness and caveats 

How robust is prediction of the dip and what are uncertainties? 
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Fig. 3. Spectra of Akeno-AGASA and HiRes (left panel) and with energy calibrated by dip (right panel). 



The shape of the calculated proton dip is stable relative to many phenomena included in 
calculations, namely, discreteness in the source distribution, mode of propagation (diffusive and recti- 
linear), inhomogeneities in source distribution on the scale I < 100 Mpc, local source overdensity and 
deficit, acceleration E max and fluctuations in energy losses. However, it is modified by the assumption 
of source evolution [B] and by presence of UHE nuclei in primary flux [§| (see also j2H] ). In Fig |4~1 
the dips for helium and iron nuclei are presented in comparison with the proton dip. One can see 



o 
CO 



CO 
o 



T3 
O 



10" 



10 



red shift 




\v He 










1 ,-■ 








— ' ' 3 


y =2.7 

'g 












4 


\ 2 


Q 17 10 18 


10 19 


10 20 10 



10 u 



o 
o 
CO 



CO 
o 



T3 
O 



10" 1 



red shift 




Y e 1 






2 \\ \ 


y=2.7 

'g 


\|4 \ 3 



10' 



10" 



E, eV 



10' 3 
E, eV 



10T 



10' 



Fig. 4. Modification factors for helium and iron nuclei in comparison with that for protons. Proton 
modification factors are given by curves 1 and 2. Nuclei modification factors are given by curves 3 
(adiabatic and pair production energy losses) and by curves 4 (with photodisintegration included). 

that presence of 10 - 20 % of nuclei in the primary flux upsets the good agreement of proton dip with 
observations. Below we shall study this problem in a more detailed way. 

Consider some acceleration mechanism operating in a gas with the mixed space density of 
hydrogen, n#, and nuclei with atomic number A , ua- We assume the power-law generation spectrum 



Qa(p) = K A n A p~ 



7s 



(2) 
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where 7 5 = 2.7 and minimum momentum is p™ in for nuclei A. We assume also that the total number 
of particles accelerated per unit time Q 1 ^ = Kn^, where K is independent of A and of charge number 
Z . Then one obtains for the ratio of generation rates of A nuclei and protons 



Qa{p)/QM 



{n A /n H )A^- 1 if pf n = v min T min Am N 
[n A /n H )Z^- 1 if pf* = R min Z 



(3) 



where T is Lorentz factor and R = p/Z is rigidity. We shall refer to the upper case in Eq. (jHJ) as to 
T-acceleration (it includes the shock acceleration) and to the lower case as to rigidity acceleration. 
For rigidity acceleration in single-ionized gas, Z c r = 1, one has 



Qa(p)/Q p (p) = n A /n H 



(4) 



This ionization condition is important only near threshold of acceleration p > p™ m , at higher energies 
the gas can be fully ionized 

The modification factor for the mixed composition is given by 



r] p (E) + \t]a{E) 
1 + A 



(5) 



where A = Q\ nm (p) / 'Qp nm (p)- The strongest distortion of proton modification factor rj p {E) is given 
by helium nuclei for which haI^u = 0.06, corresponding to mass ratio Y p = 0.24. Fig ED shows that 
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Fig. 5. Modification factors for the mixed composition of protons and helium nuclei in comparison with 
AGASA data. The left panel corresponds to mixing parameter A = 0.1, and the right panel to A = 0.195. 



rigidity acceleration in single ionized helium gas with mixing parameter A = 0.06 agrees well with the 
data, and in fully ionized helium gas with A = 0.195 (lower case in Eq.EJ) agrees worse. We shall remind 
that the first and the second ionization potential for helium is very high, 24.6 and 54.4 eV, respectively. 
In extreme case of neutral helium, the mixed modification factor r](E) w r) p (E). T-acceleration, and 
shock acceleration in particular, results in bad agreement of the p+He dip with observations. One 
should keep in mind the approximate character of our estimates, especially concerning the assumption 
about power-law spectrum down to p m i n . 

If agreement of the proton dip with observations is not accidental (the probability of this is 
small according to small x 2 /d.o.f.), Fig l5~1 should be interpreted as indication to possible acceleration 



G 



mechanism. It cannot be shock acceleration, which results in large A oc A 19 " 1 . The large 7 9 demands 
large -Emm in acceleration mechanism to avoid too large cosmic ray luminosity of a source. Correlation 
with BL Lacs indicates the jet acceleration (see Section 

Another option is given by photodisintegration of UHE nuclei by radiation. The relevant 
calculations have been recently performed in It was demonstrated that for the magnetized 

sources the observed flux above 1 x 10 19 eV becomes pure proton one due to photodisintegration by 
CMB radiation. We can add here that nuclei can be destroyed at lower energies by infra-red radiation 
inside a source. An example of such model can be given by acceleration in the inner jet, where helium 
nuclei can be photodisintegrated by infra-red emission from accretion disc. 

3. Transition from extragalactic to galactic cosmic rays 

We will follow this transition going from high to low energies. The critical energy E c ~ 
1 x 10 18 eV is given by the low-energy steepening of the dip (see Fig^J), below which the more steep 
galactic spectrum becomes dominant. The data presented in Fig ED confirm this expectation. 

In the left panel the modification factor is shown in comparison with AGASA data (see also 
Fig ED for HiRes). At E < E c the experimental modification factor exceeds 1, while by definition (see 
Eq. ^) f](E) < 1. It signals about appearance of a new component at E < 1 X 10 18 eV. 

In the right panel the spectrum for rectilinear propagation from the sources with different 
separation d and with 7 5 = 2.7 is compared with AGASA data. One can see that at E < 1 x 10 18 eV 
the calculated extragalactic spectrum becomes less than that observed. In the lower panel the similar 
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Fig. 6. Appearance of transition energy 
E c 1 x 10 18 eV in the modification fac- 
tor compared with AGASA data (left panel) , 
in the spectra for rectilinear propagation for 
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for diffusive propagation (lower panel) 
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comparison is shown for diffusive propagation with different diffusion regimes at lower energies. The 
dash-dotted curve (universal spectrum) corresponds to the case when the separation between sources 
d — > 0. In all cases the transition again occurs at E c ~ I x I0 18 eV. 
What is the reason of this universality? 

3.1. E C<1 , E c and the second knee E 2 ^ n 

We study the transition, moving from high towards low energies. E c is the beginning of 
transition (or its end, if one moves from low energies). E c is determined by energy E eq = 2.3 x 10 18 eV, 
where adiabatic and pair-production energy losses become equal. The quantitative analysis of this 
connection is given in ^Hj- We shall give here the semi- quantitative explanation. 

The flattening of the spectrum occurs at energies E < E c , where E c = E cq /(1 + z e s) 2 and z e fr 
should be estimated as redshift up to which the main contribution to unmodified spectrum occurs. 
The simplified analytic estimate for 7 9 = 2.6 — 2.8 gives I + z c g = 1.5 and hence E c m 1 x I0 18 eV. 
In fact, the right and lower panels of Fig IB~1 presents the exact result of calculations of this kind. 

In experimental data the transition is searched for as a feature started at some low energy 
E2kn - the second knee. Its determination depends on experimental procedure, and all what we can 
predict is i?2kn < E c . Determined in different experiments -E^kn ~ (0.4 — 0.8) x I0 18 eV. 

The transition at the second knee appears also from consideration of propagation of cosmic 
rays in the Galaxy (see e.g. [2%]). 

Being thought of as purely galactic feature, the position of the second knee in our analysis 
appears as direct consequence of extragalactic proton energy energy losses. 

3.2. KASCADE data and transition at the ankle 

Ankle at energy E a « 1 x 10 19 eV is seen as flattening of the spectrum in data of AGASA, 
HiRes and Yakutsk detectors in Fig ED In many works [2E] it is considered as a position of transition 
from galactic to extragalactic cosmic rays (see also [3U] for general analysis of the transition). The 
KASCADE data jH3 & ve a crue f° r understanding of transition. They confirm the rigidity model, 
according to which position of a knee for nuclei with charge Z is connected with the position of the 
proton knee E p as Ez = ZE P . There are two versions of this model. One is the confinement-rigidity 
model (bending above the knee is due to insufficient confinement in galactic magnetic field), and the 
other is acceleration-rigidity model (E max is determined by rigidity). In both models the heaviest 
nuclei (iron) start to disappear at E > Ep c = 6.5 x 1 16 eV. How the gap between I x I0 17 eV and 
I x I0 19 eV is filled? 

Protons start to disappear at E > 2.5 x 10 15 eV. Where they came from at E > 1 x 10 17 eV 
to be seen in the Akeno detector? 

To ameliorate these problems, the authors of [22] shift the ankle to 3 x I0 18 eV or assume 
another galactic component, which appears after the iron knee. The second- knee model suggests a 
solution as transition (being completed) at E c « 1 x I0 18 eV. 

The ankle model needs acceleration by galactic sources up to I x I0 19 eV (at least for iron 
nuclei). The second knee model ameliorates this requirement by one order of magnitude. 

The second-knee model predicts the spectrum shape down to I x 1 18 eV with extremely good 
accuracy (x 2 /d.o.f.= 1.12 for Akeno-AGASA and x 2 /d.o.f.= 0.908 for HiRes). In the ankle model 
one has to consider this agreement as accidental, though such hypothesis has very small probability, 
determined by \ 2 cited above. 
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E, GeV E, GeV 

Fig. 7. In the left panel the proton and nuclei spectra are shown according to KASCADE data. The arrows 
labelled by p, He, C and Fe show the positions of corresponding knees, calculated as Ez = ZE P , with 
E p = 2.5 x 10 6 GeV (the proton knee). One can notice the agreement between calculated and observed 
positions of the knees. In the right panel one can see the extragalactic proton spectrum, calculated in 
AGN model (see Section and galactic iron-nuclei spectrum (Fe-model) predicted by this model. 



3.3. Comparison of the dip and ankle transitions 

Dip transition is based on rigidity model for galactic cosmic rays, according to which the proton 
knee is at E p « 2.5 x 10 15 eV and iron knee is at E Fc = ZE P ps 6.5 x 10 16 eV. At E > E Fe the galactic 
spectrum is predicted to have steepening. The physical spectra predicted in the rigidity model are 
those for different nuclei, while the sum of these spectra has no direct physical meaning. The total 
spectrum between E p and E Fe must show some oscillations due to the knees of the different nuclei, 
and the power- law total spectrum E' 1 at E p < E < E Fe is nothing more than a fit, which cannot be 
extrapolated beyond E Fc . The question whether or not continuation of E~"< spectrum at E > E Fe fits 
the flux measured at E > E c = 1 x 10 18 eV is unphysical within this model. 

The maximum acceleration energy E m . AX for galactic cosmic rays may be in this model -Em ax < 
1 x 10 18 eV for iron nuclei, which fits SN shock acceleration models. 

At E > E Fe the galactic iron flux steepens, while extragalactic flux flattens at energies below 
E c = 1 x 10 18 eV (see curves 'Fe-model" and "extr. p' in the right panel of FigED). These two curves 
thus inevitably intersect at some energy E tT , providing the transition. The energy behaviour of both 
curves and value of E tr are model dependent. 

In Fig IT1 (right panel) the transition is shown for the model described in Section The tran- 
sition is characterised by three energies E Fe « 6.5 x 10 16 eV, E tr ~ 4 x 10 17 eV and E c ~ 1 x 10 18 eV, 
and by fluxes J ga i(£ Fe ) ~ 1.8 x 1(T 17 m-W^GeV -1 and J cxtr (£ c ) « 3.2 x 1(T 21 m-W^GeV -1 . 
It is difficult to find any fine-tuning in these values. The flux at the 'end' of galactic spectrum, E Fc , 
and the flux at the starting of extragalactic part of the spectrum, E c , differ by 4 orders of magnitude, 
and thus there is no conspiracy between these two values. 

The structure of the second knee, composed by galactic iron and extragalactic protons (see 
right panel in Fig EJ, is very similar to the structures of the other knees: there is a sharp transition 
from galactic iron to extragalactic protons (similar to transitions between different galactic nuclei), 
which in both cases results in the total spectrum with a faint transition feature. 

The predicted spectrum shape above the second knee (1 x 10 18 — 4 x 10 19 eV) is model inde- 
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pendent and is confirmed by observations. The ankle at E a ~ 1 x 10 19 eV appears automatically in 
the calculations as the part of the dip. 

The ankle model is based on the assumption that transition occurs at ankle E a ~ 1 x 10 19 eV, 
where experimental data show flattening of the spectrum. An advantage of this model is given 
by flatter generation spectrum (oc E~ 2 ' 3 ), which allows to diminish the source luminosities. The 
disadvantages are connected with models for galactic cosmic rays at E < 1 x 10 19 eV. What particles 
fill the gap between the iron knee Ep e m 6.5 x 10 16 eV and the ankle E a ~ 1 x 10 19 eV? If they 
are iron nuclei, why they have the same spectrum as the sum of different nuclei at E < Ep e 7 How 
protons, which start to disappear above the proton knee E p m 2.5 x 10 15 eV re-appear again at energy 
E > 1 x 10 17 eV comprising about 10 % of the Akeno flux? Why the spectrum at 1 x 10 18 < E < 
4 x 10 19 eV has the dip feature explained so well by extragalactic protons? 

4. The AGN model 

The AGN as sources of UHECR meet the necessary requirements: (i) to provide the necessary 
energy output, (ii) to have the space density n s ~ (1 — 3) X 10~ 5 Mpc -3 necessary for the observed 
small-scale anisotropy [321 and (iii) to provide the observed correlations with BL Lacs [53| . 

Acceleration 

According to the AGN unified model BL Lacs are Fanaroff-Riley galaxies with jets in the direction 
of the observer. Thus, correlations of UHECR particles with BL Lacs imply the jet acceleration 
of the particles. Acceleration in the jets due to pinch mechanism was suggested first for tokamaks 
(where it was confirmed on the laboratory scale) and then was rescaled in [H3] to cosmic sizes. The 




Fig. 8. Pinch-neck instability in jet 

pinch mechanism of acceleration works due to pinch- neck instabilities, illustrated in FiglHU which are 
accompanied by production of strong electric field. The solution of kinetic equations [3^ gives the 
power-law spectrum of accelerated particles q(E) oc E~ la with 7 5 = 1 + \/3- which is practically 
the same as in our calculations. The maximum energy of acceleration can well exceed 10 20 eV, if to 
rescale the laboratory pinch scale, where MeV energies were observed, to the AGN size. 

Generation spectrum and source luminosity 

Inspired by the pinch acceleration mechanism, we assume ad hoc the generation spectrum in the form: 

* { g, :; % 1 1 <«> 

with E c being a free parameter and with the maximum acceleration energy E max = 1 x 10 21 eV. 
For E c ~ 1 x 10 18 eV the emissivity needed to fit the calculations to the AGASA flux is Co = 
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3.5 x 10 46 erg/yr Mpc 3 . With space density of the sources taken from small-angle clustering n s ~ 
3 x 10~ 5 Mpc -3 , the luminosity of a source is L p = 3.7 x 10 43 erg/s, which is quite low for AGN. 

With these assumptions the calculated spectra in comparison with observations are shown in Fig ED 
The AG AS A excess, if real, needs for its explanation another component, shown by dashed line in 
the AGASA panel (e.g. from superheavy dark matter [35]). 




E, eV 



Fig. 9. The AGN-model predicted spectra compared with data. 



Transition from extragalactic to galactic cosmic rays 

Assuming that galactic cosmic rays at E > 1 x 10 17 eV are dominated by iron nuclei, we calculate this 
flux subtracting the extragalactic proton flux, calculated in AGN model, from the measured Akeno 
flux. The result is shown in Fig [7] (right panel). 

5. Conclusions 

1. The dip is the spectrum feature produced by interaction of UHE extragalactic protons with CMB 
radiation. In terms of modification factor it is practically model-independent, depending weakly 
on generation index 7 9 , mode of propagation (rectilinear or diffusive propagation), discreteness and 
inhomogeneity in source distribution, local overdensity or deficit of the sources, maximum acceler- 
ation energy and fluctuations in energy losses. For non-evolutionary sources, when the number of 
free parameters is minimal (jg and total flux normalization) the agreement with observational data of 
Akeno- AGASA, HiRes, Fly's Eye and Yakutsk (see Fig[2TJ) is very good and characterised by x 2 /d.o.f. 
equal to 1.12 for Akeno- AGASA and 0.908 for HiRes. It implies very low probability of accidental 
agreement. This is the strong evidence that majority of primaries observed at 1 x 10 18 — 4 x 10 19 eV 
are extragalactic protons propagating through CMB. 

2. The dip is modified by evolution of the sources and by presence of extragalactic nuclei as the 
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primaries. The latter effect points to some acceleration mechanism or to the sources, where He-nuclei 
are photodisintegrated (see Section O • 

3. The dip has two fiattenings, at low energy E c ~ 1 x 10 18 eV and at high energy E a w 1 x 10 19 eV. 
The latter, ankle, is confirmed in most observations (see Fig|2!J). The former provides the dominance 
of more steep galactic component at E < E c . There are three evidences of the transition from extra- 
galactic to galactic cosmic rays at energy E c shown in Fig IHU The observed transition should be seen 
at energy E < E c , and this feature is observed as the second knee. The energy E c , where transition 
from galactic to extragalactic cosmic rays is completed, is directly connected with energy E eq , where 
adiabatic and pair production energy losses for extragalactic protons are equal. 

4. The transition from galactic to extragalactic cosmic rays at the ankle is discussed in subsections 
13.2.1 and 13.3.1 For the ankle transition one must assume that agreement of the dip with observations 
is accidental, while in case of the dip the ankle is an automatic feature (part of the dip). 

5. The above-listed conclusions are valid for model- independent analysis with basically two assump- 
tions: the proton spectrum is power-law at E > 1 x 10 18 eV and sources are distributed uniformly 
with arbitrary separation d. The AGN model described in Section |^] uses some specific generation 
spectrum, inspired by pinch acceleration mechanism. This model explains small-scale anisotropy and 
correlations with BL Lacs, predicts spectra of UHECR measured in Akeno-AGASA, HiRes, Fly's Eye 
and Yakutsk experiments and quantitatively describes the transition from extragalactic to galactic 
cosmic rays at the second knee. The latter is quite similar to transitions from one type of nuclei to 
another in the first knee. 
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